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- ■ Abstract 
^ ■ 

We have established the exact expression for the gravitational potential of a homoge- 
neous polar cell — an elementary pattern used in hydrodynamical simulations of gravi- 
tating discs. This formula, which is a closed-form, works for any opening angle and radial 
extension of the cell. It is valid at any point in space, i.e. in the plane of the distribu- 
tion (inside and outside) as well as off-plane, thereby generalizing the results reported by 
Durand (1953) for the circular disc. The three components of the gravitational accelera- 
tion are given. The mathematical demonstration proceeds from the incomplete version of 
Durand' s formula for the potential (based on complete elliptic integrals). We determine 
first the potential due to the circular sector (i.e. a pie-slice sheet), and then deduce that 
of the polar cell (from convenient radial scaling and subtraction). As a by-product, we 
generate an integral theorem stating that ''the angular average of the potential of any cir- 
cular sector along its tangent circle is ^ times the value at the corner". A few examples 
are presented. For numerical resolutions and cell shapes commonly used in disc simula- 
tions, we quantify the importance of curvature effects by performing a direct comparison 
between the potential of the polar cell and that of the Cartesian (i.e. rectangular) cell 
having the same mass. Edge values are found to deviate roug hly hke 2 x lO'^ x A^/256 
in relative (A^ is the number of grid points in the radial direction), while the agreement is 
typically four orders of magnitude better for values at the cell's center. We also produce 
a reliable approximation for the potential, valid in the cell's plane, inside and close to the 
cell. Its remarkable accuracy, about 5 x 10~^ x N/256 in relative, is sufficient to estimate 
the cell's self-acceleration. 



1 Introduction 

The determination of the potential and forces due to various distributions of matter is an im- 
portant part of theoretical and numerical geophysics and especially astrophysics. Applications 
concern for instance the dynamics of solid particles (satellites, asteroids, circumstellar material, 
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planetary rings) orbiting dense objects, the internal structure of self-gravitating systems like 
plan ets, rotating stars, galaxies and massive discs, mornents of inertia and geodetic problems 
(e.g. Binnev and Tremain j . 'l987'; ^Arnab oldi and Spa rke^. 'l994'; ^Subr a nd Karas , 20051 : Nelson , 



20061 : iBaruteau and Masseti . i2008i : iVogt and Letelierl . l2009 : Schul z, 20Q9) . In general, symmetry 
in shape or in mass distribution offers a simplified mathematical framework. However, this is 
not always sufficient to enable a full analytical calculus. Other assumptions are therefore some- 
times invoked, like the reduction of the number of dimensions. The circular disc as a flat, two 
dimensional system is of great interest in astrophysics. Several cases have been studied, differing 
mainly by the size of the system (finite or not), and by the surface density profile S, generally 
a function of the radius a and polar angle, the rn o st basic being probably t he homogeneous dis c 
fiDurandl . Il958t Ihass and Blitzerl . Il988l : IConwavl . l200d IPukushimal . I2OIOI : IXresaco et al boilh . 
Non-homogeneous discs have been widely considered ii i the context o f galactic dynar nics. There 
are ax ially symmetrica l discs, a nalysed for instance by Mestel ( 1963 ) for E oc 1/a, by Casertanj 
(Il983h for S (X e"". bvlSchulzl fl2009l l2012h for S oc (1 - a)"/^ (n integer), bv lHure etal (2007^ 
and iBrada and MilgromI f ll995l ) for E oc a* (with s real). Potenti a l/surfa c e deii si ty pa i rs hav e 
also been produce d for f ully inhon aogeneous discs Kalna.is f l976h : Quanl (1993); Qian ( 1992 ): 
Evans and Collett f 1993[ ) (see also Binnev and Tremainel . 1987 ). 



In models and especially in numerical simulations working in polar and cylindrical coor- 
dinates, circular discs are naturally discretized into a large collection of massive polar cells, 
which individually generate in space an attractive force. This enables to introduce any kind of 
inhomogeneity from individual homogeneous cells. To our knowledge, the general formula for 
the gravitational potential of a polar cell — leading to these forces — is not known yet, p roba- 
bly because it is not really easy to derive from direct integral calculus. iHure et all ( l2009l ) have 
addressed the problem in the special case where the potential and field is required along the 
symmetry axis of the homogeneous cell, which greatly simplifies the mathematical derivation. 
Here, we generalize this result and report the exact expressions valid inside the distribution 
as well as outside, including the complicated off-plane case. These expressions work for any 
cell shape (radially or angularly elongated). These new results are therefore appropriate to 
estimate correctly the potential and internal forces of self-gravitating discs as continuous (i.e. 
non-particle) systems. 

The paper is organized as follows. In Section [21 we recall t he formu l a for the potential and 
acceleration of a fiat and homogeneous circular disc derived by iDurandl ( 1l953l ) in the context of 
electrostatics. In Section [3l we write the integral expression for the potential of a circular sector 
(and polar cell) and explain our strategy to perform this integration analytically from Durand's 
formula. In Section HJ we deduce the generic formula enabling to treat the circular sector and 
polar cell as well. We give in Section |5] the expressions enabling to construct the 3 components 
of the gravitational acceleration. In Section [6l we address a few remarks about technical 
aspects (computation in the cell's plane, potential in the polar axis, long-range behavior). The 
formula for the potential is graphically illustrated in Section [71 In Section [SI we generate a 
theorem about the angular average of the potential in the cell's plane, deduced from periodicity. 
Section [91 is devoted to a comparison of the potential of the polar cell with its Cartesian 
analog, in conditions typical of current numerical simulation of hydrodynamical discs where 
polar cells have almost rectangular shapes. This helps to decide when and where curvature 
effects can be neglected. Finally, in Section [TDl we establish a very good approximation for the 
po tential, by using th e expansion of the elliptic integral of the first kind around the singularity 
bv lVan de Veil f ll969l ). A few concluding remarks follow. Several appendices contain definitions 
and demonstrations. 
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2 Potential of a flat, circular disc 



The derivation of the expression for the potenti al of a homoge neous and flat, circular disc, 
with radius a as shown in Fig. [1^ is found in iDurandl (119531) in his rema r kable textbook 
devoted to Electrostatics (for the astrophysical context, see iLass and Blitzerl . Il983l : IConwavl . 
2000 : Tresaco et al . 201ll ). Up to a factor — GS where G is the constant of gravitation and S 
is the surface density of the disc, the expression reads, retaining Durand's notation: 



disc / 



r: a) 



a + K 



where R and Z are the cylindrical coordinates of the observation point P{r) (i.e. r = Rcr + 
Zez), E, K and 11 are the complete elliptic integrals of the first, second and third kinds 
respectively (see the Appendix |X] for a definition). 



where z is the altitude of the disc, 



c 



(2) 



is the modulus (0 < A; < 1), 



5=^(a + Ry + C, 



k 



S 

2\/aR 



m 



{a + R) 

is the characteristic or parameter (0 < < m < 1), and 



(3) 
(4) 

(5) 



1, if i? < a. 



1 

2' 



if R 



0, if R> a. 



(6) 



which is a Heaviside-type function, i.e. e'+H(i? — a) = 1. he potential is finite and continuous 
everywhere, including at the edge (i.e. at i? = a and Z = z). I t is possible to express the 



Il-function in terms of Heuman's lambda function (jConwavi . l2000l ). 

We can easily calculate the potential everywhere in space of an annulus, simply by consid- 
ering Eq.([T]) with two different values of a. For instance, if denotes the inner edge and Oout 
is the outer edge of the annulus, then we have: 



annulus / 



r, O-in, C^outJ 



disc/ 



(7) 



Equation ([T]) is also the starting point to determine the potential of any axially symmetrical 
system with non-zero thickness, by performing an integration along the 2;-direction. This is 
usually achieved through a numerical process. Finally, from Eq.(IT]), we can get the gravitational 
acceleration g'^^^'^ = —'Vtp'^^^'^ by differentiation (see appendix[B]), which is an important quantity 
to study the dynamics of test particles. Note however that g'^^^'^ logarithmically diverges at the 
edge (like for any homogeneous, bidimensional distribution). 
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Figure 1: Configuration for a circular disc witli radius a {top panel,) and for a circular sector 
with radius a, bounding angles 6'^ and (negative here), and opening angle A6'' {bottom panel). 

3 Relaxing the axial symmetry: strategy of calculus 

We now consider a homogeneous circular sector as the material domain bounded by two polar 
angled O'^ and ^2 ^ ^'i? the polar axis i? = 0, and the circular arc with radius a. The corre- 
sponding configuration is shown in Fig. [T]d. The opening angle is A^' = O^ — O'l, and the mass 
(per unit surface density) is: 

M = ^a^Ae'. (8) 

In the following, we assume, without loss of generality, that 6[ and 62 are in the range [— vr, tt]. 
This enables us to describe any circular sector, ranging from the line (the asymptotic case where 
6*2 — )■ 9[) to the circular disc (the case with 62 = —9[ = 71). To get the a ssociated p otential in 
space, it is necessary to go back to the integral expression which i^ (e.g. Zhu, 2005 ): 



V'-^*-(r; a, ^,,^2) = £ [F(/3i, k) - F(/32, k)] da, (9) 

where the amplitude /3 of the incomplete elliptic integral of the first kind F (see the Appendix 
|A]for a definition) is defined by: 

T,-{e' -e) = 2/3, (10) 

and 9 denotes the polar angle of the observation point P (see note [1]). We then have 2/3i = 
TT — {6[ — 6) and 2/^2 = 'n — {62 — 0). Since /32 < /^i, the above integral is always positive. With 



"'^As usual, polar angles are measured positively from e^- 

^To lighten the notation, we use here the same variable a for the integration variable and for the integral 
upper bound. 
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P(R,9,Z) 



(a) floating sector 





Figure 2: A "floating" circular sector {top panel) is a particular circular sector with radius a and 
bounding angles 9[ = 9' and 62 = 9 + tt (this second angle depends on the observation point, and 
so the floating sector is not very physical). Its opening angle is A^' = 6 + it — 6'. Conflguration 
for a polar cell {bottom panel) with inner radius ain, outer radius Oout, and bounding angles 9[ 
and 




the above deflnition, we see that the quantity: 

^kF{f3,k)da = ^^°''\f;a,P) (11) 
R 

represents the potential of a circular sector deflned by the two bounding angles 6[ = 9' and 
9'2 = 9 + n (i.e. /3 = 0). Such an object is not very physical since the second bounding angle 
depends on the polar angle 9 of the observation point P (see below for 9 = 9'); it is introduced 
for convenience as it plays a central role in the calculations throughout. We call such a circular 
sector a floating circular sector. It is displayed in Fig. [2^ (in the case where /3 > |). The 
potential of a circular sector is then given by the difference of the potential of two floating 
sectors, namely: 

r^^'°\r; a, (3„ P2) = V^'°^*(r; a, P,) - i^^^^^r; a, P^). (12) 

We note that for 6*2 = —9'^ = vr = (i.e. the circular disc), we have 2/3i = 2-n — 9 and 2/^2 = 9. 
In these conditions, we get: 

F(/3i, k) - F(/32, k) = F{n + ^, A;) - F(^, k) (13) 

= F{7r,k) + F{^-,k)-F{^-,k) 
= 2K{k) 



5 



and so Eq.(j9]) takes the well-known form: 



^'^'\r-a) = 2l J^kK{k)da, (14) 

Jo \ R 

whose closed-form is therefore Eq.([T]). In the special case where /? = f (i.e. 6 = 6'), the floating 
sector is exactly a half-disc, and the observation point P is somewhere along the diameter. As 
a consequence, ip^"^^ is half the potential of the circular disc, i.e. ^ip'^^^'^ {r; a). For amplitudes 
/3 7^ ^, the determination of Eq. ffTT]) by analytical means is not trivial. We can distinguish two 
cases: 

• within the plane of the distribution where ( = 0, we have k = m. This is the most simple 
case: the result can be found either by a change of integration variable {da — )■ dm), or by 
using a change of mod ulus of the . F-fun ction (i.e. m— )-a/i?<lorm— )■ R/a < 1). This 



is the method used in iHure et all (120091 ) to determine the potential along the symmetry 
axis of the polar cell. 

• when the observation point is off the plane of the sheet (i.e. C 7^ 0); things are more 
tricky (the goal of the paper). We have not been able to solve this general problem by 
direct integration but by exploring the incomplete version of Durand's formula. It turns 
out that this approach was successful. 

To get the incomplete version of Durand's formula, we simply change the complete elliptic 
integrals into incomplete ones. For each term in Eq.([T]), the correspondence is then: 



6E{k) -> 5E{(5,k), 

^K(fc) ^ «!^F(/3,fc), (15) 
f^^^(m^A:) ^ ^^n(/3,m^fc), 

where E and 11 are the incomplete elliptic integrals of the second and third kinds respectively 
(see Appendix |X] for a definition). We have then calculated the partial derivatives of each 
term with respect to the radius a, aiming at the generation of the right integrand ^^J^kF{P, k) 
appearing in Eq.( lTT|) . This process requires the calculation of the partial derivatives of the 
incomplete elliptic integrals with respect to k and/or m (see Appendix [C]), as well as the 
derivatives of the modulus k and parameter m with respect to a, namely: 

and 

— d Tfi = (17) 

m " 2a{a + R)' ^ ' 

The derivatives of the three terms in the right-hand-side of Eq. (fT5|l are fully expanded in 
Appendix ini By summing the resulting formulae, we find that i) the incomplete elliptic integrals 
of the second and the third kinds disappear, and ii) the expected integrand is indeed present. 
We actually obtain the following relation: 



da 



SE{P, k) + ^--^F{/3, k) + ^^n(/3, m^ k) 
a + K 

-kF{/3,k) + — — sm/3cos/3- . . 2 o ■ 

K 2a 1 — sm p 



;i8) 
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and we recognize the expected integrand ^y^kF{l3, k). The point is that the second asymmetric 
term can be written in the form of a partial derivative with respect to a too. This was a priori 
not guaranteed, but there is indeed a closed-form. The demonstration is given in Appendix [El 
For any angle /3 < |, we have the general relationship: 



m?5 sin /3 cos — sin^ 13 



2a 



[1 — m? sin^ (3) 



da { -R sin 2/3 asinh 



+Carctan 



a + i? cos 2/3 
VC^ + R'^ sin^ 2/3 
C(a + i?cos 2/3) 



(19) 



R sin 2/3 - 4a/? sin^ /3 



and so we deduce the indefinite integral: 



' ^A;F(/3, k)da =6E{P, k) + ^F(/3, k) + ^^n(/3, m^ fc) 
K a + K 



(20) 



— sin 2/3 asinh 

— Carctan 



a + -Rcos 2/3 

VC^ + sin^ 2/3 
C(a + i?cos2/3)' 



where 



pp/ 



sin 2/3 PP' 



52 -4a/?sin2/3, 



(21) 



is the distance from point P(-R, 6, Z) to point P'(a, 6', z). To be rigorous, we should include in 
the right-hand-side of Eq. fl^U]) any function of the three variables /3, and i?, but this is not 
necessary here since we work with a definite integral in the following. With and a as the 
lower and upper bounds respectively, we get the expression for V^*'°^*(f) defined by Eq.([II]), and 
then ^'^e^t^ff) given by Eq.([9]) or Eq.([l2]). As ^(/3,0) = F(/3,0) = n(/3,0,0) = /3, we deduce 
that the elliptic integrals cancel each other at the lower bound. We finally get, for /3 < |: 



^^°^*(f;a,/3) 



5E(/3, k) + ^F(/3, k) + ^$n(/3, m^ A;) 
a + it 

+RH{P) + CT{f3), 



(22) 



where we have defined for convenience: 



r(/3) 



-arctan 



C(a + i?cos 2/3) 



/?sin2;g PP' 
for the "trigonometric" term, and 



+ arctan 



c 



cot an 2/3 



(23) 



sin 2/3 asinh 

V vC 

— asinh 



a + R cos 2/3 
cos 2/3 



(24) 



+ sin^ 2/3 



for the "hyperbolic" term. Note that both H and T depend not only on the amplitude /3 but 
also on m and k through a/R and C,. Equation fl22l) is exact. It is valid in all of space, in 
the plane of the distribution (if ^ = 0) as well as off the plane. It can obviously be put in 



7 



various equivalent forms, in particular by combining the intrinsic functions having two different 
arguments. As expected, it is an even function of ( (the potential is the same above and below 
the plane of the distribution). 

Because of the inverse trigonometric function in Eq. ffTOj) . Eq. fl22|) must not be used for 
amplitudes (3 larger than | (unless T must include a constant). Such amplitudes correspond 
to a floating sector larger than half-a-disc, or A6' > tt. This configuration is precisely the one 
displayed in Fig. [2^. To treat cases with /3 > f, we replace this large floating sector by a circular 
disc minus a small floating sector. The small floating sector to be removed has the following 
properties: i) same radius a, ii) opening angle 2n — A6' < tt, and iii) two bounding angles 
e" = 26-6' >e andO' = 6 + 71 > 6". The corresponding amphtudes are /3" = f - ^ < f 
and /3 = 0, respectively. For /? > f , the potential of the floating sector is then given by: 

^^°^\f; a, /3 > |) = ^^'""'{f; a) - ^^°^\f; a, (5" < ^). (25) 

A typical algorithm to determine ip^"'^^ whatever /3 by using uniquely Eq.f E^ is then the 
following: 



let a^6' -6- 

let/3 ^ 

if /3 < f then 

let result ^ ^float^- 

else 

let/3' ^ 

let result ^ 7/)fl°at^- ^/^. 
let result ^ 2ilj^°''^(f] a, |) - result] 
endif 



Figure |3] displays ■^/j^o^.t^-. computed in the (-R/a, /3)-plane from Eq.f l22|) . for R/a E [0, 2] 
and (3 e [0,7r]. When /3 > |, Eq.f E^ is used instead. Two cases are considered : C = (i.e. in 
the plane of the circular sector), as well as an off-plane case, where C, = a/2. 



(26) 



4 Potential of the circular sector and polar cell 

The potential of a circular sector is found from Eqs. fl2^ and fl2^ by considering two amplitudes 
/3i and P2, according to Eq. (|T2|) . Similarly, the gravitational potential of a polar cell is found 
from Eq. (|T2i) by considering two circular sectors having the same bounding angles but two 
different radii. A typical polar cell is depicted in Fig. |2]d. Let Oin be the inner radius of the 
cell, Oout the outer radius, Aa = Oout — > the radial extension, and 6[ and 62 > 6[ the two 
bounding angles. So, the potential is: 

r'\f; a;,, ao„t, /3i, /32) = a^ut, /3i, /32) - r'''°'{r; ai„, /3i, fS^), (27) 

where, for each circular sector, the quantities 6, k and m must be computed as: 



(ain + ry + e 



2y/a~R 
5 ' 



(2J 



for a 
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0.5 1 1.5 2 

normalized radius R/a 



Figure 3: Color map showing tp^°^^{f,a,P) from Eq.f E^ versus R/a and /3 (in units of vr). 
These cases are for C = {top panel) which corresponds to the plane of the floating circular 
sector, and for ( = a/2 (off-plane; bottom panel). The color code is the same. A few contour 
levels are superimposed 
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and 



k = 2v^^^ (29) 

/ 

3 

while ( is the same for the two sectors. 



^ = llf^' for a = ao,t, 



5 Accelerations 

The gravitational accelerations —S/ip'^'^'^^"'^ and —Vip'^^^^ due to the circular sector and polar cell 
respectively involve the partial derivatives with respect to the three variables R, 9 and Z. We 
have 

Wr'^°\r; a, = VV>fl°^*(f; a, ft) - V/°'^*(f; a, /^s), (30) 

and 

VV^^^"(r; ai„, Cout, A, P2) = V i^^'^^'^'^r- a^ut, A, (32) - V^'^'^'ir; a,^, (3,, /^a), (31) 

where — V'?/'^°'^*(r; a, /3) is the acceleration of the floating circular sector with bounding angles 
and 6' = 7T + 6 — 2/3. The 3 components of V'?/'''°'^* (r*; a, /3) are derived in Appendices [0 [G] 
and [HI As for the potential (see Eq. (l25l) ). the case with /? > f must be treated by considering 
the circular disc, i.e. 

V^^°'^*(f; a,(3>^) = V^'^'^^if; a) - V^^°^\f; a, (3" < ^), (32) 
where /3" = vr - /3 < f , and V^'^^"'=(f; a) = -^discj^- ^^^^ Appendix E]) . 



6 Important remarks and special cases 

Case with /3 = |. As mentioned above, this case corresponds to the value of the potential 
at the edge of a half-disc (i.e. along the axis defined by 6* = 9'). For that amplitude, all 
the incomplete elliptic integrals take their complete value. The inverse hyperbolic functions 
vanish, while the two inverse trigonometric functions may still contribute through an argument 
eventually infinite. We then need to consider three cases, depending on the position of the 
observation point relative to the tangent circle, which is the circle with radius a and tangent to 
the circular sector. Th e same kind of distinction exists under axial symmetry fiDurandl . Il953 : 



Lass and Blitzerl . 119831 ). We have: 



if R 



then m = 1 and a + i? cos 2/3 = 2i?cos^/3. The argument of the first inverse 



tangent function remains finite. We find: 



TT 



T(/3) = — Id arctan (cotan l3) + — 



(33) 



and so 



hm T(/3) = -|CL 
2 



(34) 



if i? > a, then {a — R)/ sin 2/3 < 0. The two inverse trigonometric functions cancel each 
other, and so we get: 

lim T(/3) = 0. (35) 
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ii R < a, then {a — R)/ sin 2/5 > 0. The two inverse trigonometric functions bring the 
same contribution and we find: 



hm T(/3) = -ICIvr. (36) 

We then conclude that: 

I 1, if R < a, 

(37) 



1, 


if i? < 


a. 


1 

2' 


if R = 


a. 


0, 


if R> 


a. 



We see that this result can be written in the compact form: 

lim T(/3) = -ICItt X e', (38) 

where e' is a step-function already defined in Section [2j We then fully recover Durand's formula. 

Axial symmetry. Given the above analysis, we see that Durand's formula and Eq. (fT2|) are 
indeed equivalent. Actually, axial symmetry implies = 271, which is achieved for instance 
with 6[ = 9 and 62 = 9 + 271, that is for 2/3i = vr and 2/^2 = —71. In particular, the term 
— |C|vr X e' in Eq.([T]) comes from the two trigonometric functions, in the limit /3 — We 
have: 

^'^'-(r; a) = V^-^*-(f; a, |) - ^^'^'"^{f- a, -^) (39) 
= 2^fl-Hf;a,^), 

which is nothing but the potential of two half-discs, hence the factor 2 in Durand's formula. 

Calculations in the plane. For ( = 0, we have k = m and 6 = {a + R). The two functions 
n and T do not contribute. Equation (ITT]) then becomes: 

^fl°^t(r = ReR-a,P) = {a + R)E{/3,m) + {a-R)F{/3,m) (40) 

o • o« ■ u « + ^cos2/3 . COS2/3 
— Rsm2p asmh — — — — asmh 



V /2|sin2/3| \sin2/3\J ' 

Polar axis. We see from Eq.f l40p that for R = and ( = 0, S = a, and k = m = 0. Only the 
incomplete integral E and F contribute, and they are identical to the amplitude. The potential 
of the fioating sector (up to a factor — GS) is then given by: 

?A'i°'^*(0;a,/3) = 2a/3, (41) 

which is nothing but the length of the arc. So, the potential of the circular sector at its corner 



is: 

For ( 0, we have 
and then 

which leads to: 



V'^^^*°'^(0; a, A, (32) = 2a(/3i - /^a) = a AO'. (42) 
T(/3) = -2/3|C|, (43) 



^fl-t(f = Zez;a,(3) = 2/3 (^^^^Te - iClj , (44) 
^sector^^-.^ Zez;a,/3i,/32) = [V^FTC - |C|) A^'. (45) 
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Long-range behavior. At large distance from the system, the modulus A; — )■ 0, meaning 
that F{/3, k) ^ P to second-order in k. At the same order, we then have 

V^^°^*(f;a,/3) ^ / —Pda. (46) 

./n d 



With S ^ r ^ ya? + z^, r being the spherical radius, we have 

^^°^*(f; a, /3) ^ - r '^da = -(3a^. (47) 
r Jo (3 r 

By difference, we get from Eq. fll2p : 

M 

r'^'^'^'^if; a, /32) ^ -(/3i - P2) = = -, (48) 

which is the expected result: far enough from the system, the potential is proportional to the 
mass and inversely proportional to the relative distance. 



7 A few numerical examples 

The way to compute the potential depends mainly on the polar angle 6 relative to 6[ and 62. 
We recommend to proceed as follows: 

• if ^2 > ^'1 — ^5 ^he two amplitudes /3i and P2 are less than ^. The potential is then 
computed directly from Eq. f|T2|) where each term is determined from Eq.f l22|) : 

• if 6[ < < 9, we can proceed as above by considering a new observation point P"(F') 
mirror symmetric of P with respect to the plane of the distribution for which the potential 
is the same. The altitude of this point is 2z—Z (the radius and polar angle are unchanged). 
The new bounding angles are : 9'1 = —O^ and 62 = —6[- The potential is then 

r^^'^^f; /32) = r'^'^'^'^if"; f3[, /3^), (49) 

where 2/3' = n - {9" -6). 

• i{e'2>9 > 9i, then /32 < f < In this case, tp^°^\f] a, P2) is determined from Eq.(|22D 
and ■?/;fl°^t(f; a, (3i) is found from Eq. fl2^ . We finally have (see also Eq. fl2^ ): 

^sector ^ ^disc^^. _ ^Aoat^^. ^z^, _ ^Aoat^^. ^^^)_ ^gg^ 

where (3[ = n — (3i. 

Figure H] shows the potential of a circular sector in two [R, 6')-planes. The first panel is for 
C = (i.e. in the plane of the cell), and the second panel is for ( = a/2. The opening angle 
of the sector is A9' = |, i.e. 9[ = —^2 — ~j- Finally, the potential of a polar cell computed 
from Eq. fl27|) is shown in Fig. Elfor C = and for ( = a/2 (off-plane). The opening angle is the 
same, the radius at the outer edge is Oout = foin and Om = 1. 
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Figure 4: Color map showing the potential 7^^^'=*°'' of the circular sector (the boundary is shown 
in white) with opening angle A6' = |: in the plane {top panel) and off-plane with ( = a/2 
{bottom panel). The color code is the same. A few contour levels are superimposed. 
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Figure 6: The odd function H{(3) for a = i? and C = versus /3 (in units of n). 
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8 A theorem 

In the plane of the distribution, Eq. fl22]) takes a particularly simple form since m = k (i.e. 
C = 0). In addition, for R = a we have m = k = 1. All along this tangent circle, we have: 

^fl°^t(r = aeR;a,/3) = 2aE{(3,l) + RH{(3), (51) 

where 

HiB) = — sin 2/3 | asinh — ^ — asinh ^ | , (52) 

|sin2/3| |sin2/3|y' ^ ' 

= -sin2/3 X In ( 1 + ^ 



cos/3| 



The function H{(3) is plotted versus /3 in Fig. [6]for a = i? and C = 0. It is an odd function of 
the angle 2/3, and is 2-7r-periodic in this same angle. As a consequence, when integrated over 
any interval [0,7r], we always get: 

H{l3)d(3 = 0. (53) 





In fact, this periodicity does not depend on a and R, and so Eq. fl53l) is true even for m 7^ 1. 
We can easily integrate over one period each term in the right-hand-side of Eq. f|T2l) . We take 
the interval [9'i — tt, 9[ + n], which means 2/3i G [0, vr]. As we have 

'sin/3 if/3G[-f,f], 
E{f3, 1) = { (54) 
2-sin/3 if/3 G [? 



L2 ' 2 J' 



we find: 



iP^°^\f=aeR;a,/3i)d9 = 2 / ip^°^\f = acR] a, (3i)d/3i (55) 

Jo 

= Ana. 



Regarding the second term, the integration which is to be performed over the same interval in 
9, gives: 

^fl°at(f = aen; a, /32)d9 = 4a(7r - A^')- (56) 
Consequently, the potential of the circular sector integrated over the whole tangent circle is 

^sector aen; a, f3i, /32)d9 = 4aA9'. (57) 

27r 

We can therefore postulate that the angular average of the potential along the tangent circle of 
any circular sector with radius a and opening angle A9' is (up to a factor — GS): 

— / ij'^^'°'{9)d9 = 4af, (58) 
27r Jstt 

where 

/ = ^ e [0, 1] (59) 

/TT 
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Figure 7: The potentials il^^"^'^ (^r = ae^] a, and tp^'^^^{j = ae^] a, (32) as well as their difference 
^sector ygj^gijg along the tangent circle {top panel). Here, the parameters are a = 1 and 
6*2 = f = —9[, and so A9' = | and f = \- According to Eq.f l57|) . the area (shown in grey) 
under this last curve is a. The potential versus the polar angle 6 along the circle tangent 

to the circular sector, for different opening angles A^' labelled on the curves, and a = 1 (bottom 
panel). The area under each curve is iaAO' according to Eq. flS7|) . 



is the opening factor of the circular sector. This relation can be rewritten into different forms. 
In particular, by using the value at i? = (see Sect. [6]), we also have: 

— / ij'^^'°\e)de = -V''''*°'(0). (60) 

meaning that the angular average of this potential is ^ times the value at the corner. For the 
circular disc, we have A^' = 2tt and so 



1 

2^ 



^lj^''''{e)de = 4a. (61) 



27r 



Figure [7^ shows 7/;^°^* versus 9, as well as the potential of the corresponding circular sector 
^sector_ Figure [7)d shows ip^'^^^"'^ versus 9 over one period for various opening angles A^'. 

It is probably possible to derive other theorems by considering only the range [9'^., 9'<^ (instead 
of a full period), or integration along the radius. 
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9 Curvature effects and high resolutions 



When the opening angle AO' of the polar cell tends to zero, curvature effects are expected to 
gradually vanish. The cell becomes more and more rectangular, i.e. Cartesian. It is interesting 
to analyze how potential values vary in such a case. In particular, it can be interesting to be 
able to decide whether or not the formula for the Cartesian (i.e. rectangular), homogeneous 
cell can be used instead of the formula for the polar cell, for a given accuracy, and to quantify 
the deviation between the two. This is what we have done by considering practical cases. In 
current numerical simulations of self-gravitating discs, the discretisation consists in a grid with 
N X L cells (A^ is for the radi al direct i on, and L is for the angular direc t ion), typically a few 



hundred in each direction (e.g. iNelsonl . l2006l : iBaruteau and Massetl . 120081 : iPierens et all . l201ll ) 



The angular resolution j is often smaller than the radial resolution in order to produce, at 



least at some radius, polar cells with almost square shap^. This is achieved when 



aAe' ~ Aa, 



(62) 



where A9' is the opening angle (see Sec. [3]), and Aa is the radial extension (see SecS]). Since, 
by nature, cells have not the same size in polar coordinates, while the shapes stay homothetic, 
this conditions cannot be met everywhere. For a regular radial discretisation, this condition 
translates into: 

(63) 



6A^, while at twice this distance, 



L N 

and we see that the polar cell is close to the square cell L 
this is achieved for L = 3N. 

Like a polar cell, a Cartesian cell is defined by four corners with coordinates {x[,y[), (xg, y'l 



and {x'i,y2) and its mass is M = (x'a — x[){y2 — y[), per unit s urface density. The 
expression for the potential of such a sheet, denoted ip'^^'^^', is known (e.g. Durandl . 1953 ): it is 
recalled in the Appendix [H The potential of a rectangular sheet with — x[ = y2 — y'l = is 
shown in Fig. [H^. We have then compared the potential of the polar cell ^Z'^®" and the potential 
of the Cartesian cell ip'^^^^- in the case where the two systems have the same mass, for various 
numerical resolutions 1/A^ and 1/L. To complete the above example, we show in Fig. [8]d the 
potential of the polar cell having the same mass and same radial extension as the rectangular 
sheet. We have compared the two potentials numerically. For this purpose, the setup is the 
following: 



1, 



-'out 
^2 = 



1 + 

2tt_ 
' 2L 



N 



(64) 



for the polar cell, and 



X2 

y'2 



•'out 



M 



(65) 



2(4 -^'i) 



-y'l 



for the rectangular sheet. The value of y'-^ ensures that the mass of the two systems remains 
the same whatever the parameters N and L. The relative difference 

^cell 

A = 1 - -— (66) 



rect. 



^This preserves a certain isotropy of the precision of the differentiation schemes of fluid equations. 
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Figure 8: Color map showing the potential ip^'^'^^- of a rectangular sheet in the (x, y)-plane 
computed from Eq. fll36p {top panel). The four corners have Cartesian coordinates 
(|, 1/^), (1, —y'l) and (|, —y'l), where y[ is such that the mass of the rectangular cell is the same 
as the mass of the polar cell, according to Eq.fl65b). The boundary of the Cartesian cell is shown 
as a white line. A few contour levels are superimposed. Same legend but for the potential ip'^^^^ 
of a polar cell computed from Eq. fl271) [bottom panel). The parameters are ( = 0, Oin = 1, 
ctout = I and A9' = n/6, i.e. N = 2 and L = 6 according to Eq. flMl) . The boundary of the 
Cartesian cell with same mass is recalled (dotted line). 
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X 

Figure 9: Color map showing the relative deviation A defined by Eq. (l66|) between the potential 
of the polar cell and the potential of the Cartesian cell with same mass. The parameters are 
= 1, iV = 256 and L = QN (see also Figs. [S^ andlHb)- A few contour levels are superimposed. 

is shown in Fig. [9] for = 256 and L = 6N inside the cell and close neighbourhood. As 
expected, the largest deviations are observed at the four corners and along the boundaries, in 
excess or in default. The maximum value is A^ax ~ 2 x 10^'^ here (in absolute). The deviations 
remain globally small : i) along the mean arc R = |(ain + dout), and ii) along the symmetry 
axis 9 = ^{9'i + ^2)- This means that the radial and angular accelerations and forces for the 
rectangular sheet and for the polar cell are close to each other there. The deviation at the center 
of the cell Ac is particularly small, with Ac ~ 10^^. 

Figure [TO] shows Amax and Ac for running from 2^ = 16 to 2^^ = 4096 and L = 6N still. 
Both log |Ajnax| and log |Ac| vary as a power-law of A^. There is about 4 orders of magnitude 
between the deviation at the corners and the deviation at the center. We have roughly: 

, ^ /256\ 

|A^ax|~2x 10-3 X i—j 

for boundary values, and 

s /256\^-^ , , 

|Ac| ^4x 10-« X ( — j (68) 

for values at the cell center. We have considered different ratios L/N, in particular L = 3N and 
L = 12N which corresponds to polar cells slightly elongated in the angular and in the radial 
directions, respectively. We have noted that all the results are globally similar. 

10 Approximation in the vicinity of the cell 

As a matter of fact, we can find a good approximation for the potential in the vicinity of the 
polar cell. This may actually be interesting not to manipulate incomplete elliptic integrals and 
to work with a less complicated expresssion. The vicinity of the polar cell means /3-amplitudes 
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Figure 10: Relative deviation A between ip^'^^^ and ^|J^'^'^^■ as a function of the resolution param- 
eters N and L = QN: at the center of the polar cell {white square) and the maximum value 
{filled circles) taken over the whole potential map (see also Figl9]for = 256). 



close to I and modul i m close to unity . In such a case, we can use the following expansion 
based on the work by I Van de Veil (Il969l ): 



'E{l3,m) ^ 1, 
F{(3, m) In ■ 



4tan/3 



In 



4(2+ro) sin/ 



where m' = \/l — is the parameter complementary to m, 



(69) 



and 



W{w,l3) = ^/A{l + tu) COS2/3 + 
R-a 



(70) 
(71) 



In contrast, it does not seem possible to approximate the hyperbolic term H{f3) for m = k 
without generating noticeable errors. So we take its exact form: 



H{(3) = - sin 2/3 x In 



1 + 



W{zu,P)\ 



(72) 



The approximation for the potential of the floating sector, denoted ip^"'^^, is found from 
Eq.(l22]) or Eq.(^ using Eq.(l69]). We get: 



<°'^*(i?e,j;a,/3) = a 



2 + — ti7 In 



4(2 + tu) sin/3 
(2 + CC7) cos/S + W 



-{1 + w) sin2/31n ( 1 + ^ ) 



(73) 



For both coherence and optimization, the value of ip'^^^^, when necessary (i.e. in the case where 
/3 > |), must be replaced by its approximation at the same order, namely: 



^t"iReR;a) = 2a 



w + 2 + vj\n 



\zu\ 



A{zu + 2) 



(74) 
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The approximation for the potential of the circular sector and then for the polar celjfl are 
then found from Eqs. (175i) and (frll) by using convenient amplitudes (3i and /32, and radii and 
Qout- We have checked the accuracy of this approximation in the same conditions as above. 
The relative difference 

^cell 

is shown in Fig. [11] for = 256 and L = QN inside the cell and close neighbourhood. We note 
that the agreement is remarkable with |A| ?a 10~^. It is especially excellent again along the 
mean arc R = |(ain + ctout); which means that the angular component of the force should be 
very accurate. 

As above, we have changed the resolution (still with L = 6N). The results are displayed in 
Fig. [121 We roughly find: 

/256\ ^'^ 

lAboundaryl ~ 5 X 10"^ ^ ( ^ ) C^^) 



for values of the potential along the boundary, and 



256\^ 



|Ae| ^ 8 X 10-« X ( — ) (80) 



for values at the cell center. 



11 Concluding remarks 

In this article, we have investigated the gravitational potential of the circular sector and polar 
cell, which are two fundamental patterns present in numerical simulations of disc-like structure. 
Several new results have been established : the exact formula valid in all space and for any 
cell shape, an average theorem and an approximation for the interior potential. We have 
also quantified the curvature effects by comparing the potential of the polar cell and that of 
the Cartesian cell. All the formulae presented here are appropriate to both theoretical and 
applied studies requiring either low or high spatial resolution. Applications to electrostatics 
and capacitance problems are straightforward. 



is: 



^Inside the cell, we have Z?, < f < and so the approximation for the circular sector (i.e. for 6 E [9[, 62]) 

^'""'^°'{ReB:,a,l3i,l32) ^ a {2m In \w\ + w \n{sm j3[ sin ^2) (75) 
+Tuln {[(2 + w) cos /32 + W^(/32)] [(2 + vj) cos I3[ + Wil3[)]} 



+ (1 + n7) 



sin 2/3[ In ( 1 + f ) + sin 202 In ( 1 ^ 



where (3[ ^ + ^li^ — 6 = n — Pi. The radial acceleration is given by 



gT'°'iReR;a,Pi,P2) - --9^V'°=*"(i?efl; a, /3i, Z?,), (76) 

a 



and the angular acceleration is 



jr'°'iReR; a, /S,) = -^der''''°'iReR; a, Pi,^^). (77) 



21 




X 



Figure 11: Color map showing the relative deviation from Eq. fl78|) between the exact expression 
for the potential ip'^'^^^ and its approximation tpQ*^^^ in the same condition as for Fig. [HI A few 
contour levels are superimposed. 
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10 100 1000 
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Figure 12: Relative deviation A between ijj given by and its approximation as a function of 
the resolution parameters and L = 6N: at the center of the polar cell {white square) and 
along the boundary {filled circles) (see also Fig{TT]for N = 256). 
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A Elliptic integrals 

The incomplete elliptic integral of the first kind is 

F(0,A;) 



dx 



V 1 — gin X 

where k is the modulus and (p the amplitude. The incomplete elliptic integral of the second 
kind is ^ 

E{(p,k)= I dx\/l - F sin^x, (82) 
Jo 
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The incomplete elliptic integral of the third kind is 



n(0,m^A;) 



dx 



(1 — m? siv? x)\/l — kP' siv? x 



13) 



where m is the characteristic or parameter. The complete elliptic integrals are found by setting 
= |, namely 

(84) 



F(-,k) = K{k), 



. TT 



E(-,A;)^E(A;), 



and 



,7r 



n(^,m^ A;) = n(m2, A;). 



^5) 
^6) 



B Acceleration of the circular disc 



The two components o f the accelerat ion g^^'^'^ = —Vip'^^^^ associated with Eq.([T]) are also given 
in Durand's textbook (iDurandl . Il953l ). Up to a factor — GS, these are: 



disc 



1--)K(A;)-E(A;) 



(87) 



and 



dzi> 



disc 



27rSign(C)e' + 



2C 



R — a 

R + a 



n(m, k)-K{k) 



where Sign(C) is the sign of ( (or if C = 0). 



C Derivative of the elliptic integrals 

For the incomplete elliptic integrals of the first and second kinds, we have respectively for any 
modu lus k, complementary modulus k' = \/l — /c^ and amplitude (j) (iGradshteyn and Ryzhikl . 



20071 ) 

kk'%F{<j>, k) = E{<p, k) - k"F{<j>, k) - ^''^^'^'"'f , (89) 

V 1 — A;^ sin'' 

and 

kdkE{(l),k) = E{(l),k) - F{(P,k). (90) 
Regarding the incomplete elliptic integrals of the third kind, we have two derivatives. We have 

a,n(</., m^ k) = k -^^^^^^^^==3 , (91) 



( 1 — sin^ 0) a/1 — A;^ si 



sm^ I 



and so 



A;^ — ^ ^ , , , , , — w?) sin^ , . 



^ Jo (1 - w? sin^ 0)\/l - A;2 sin^ 

= -n(0,m2,A;)+ T — 

JQ ./l 



-^(0,m^A;) + 



A;2 sin^ (j) 
E{(j),k) A;^ sin cos I 



k'^ k'W^-k^ sin^ 
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where the last integral can be found in 
with respect to k is then given by 



Gradshteyn and Rvzhik 



(120071 ). The partial derivative 



n(0,m^fc) - 



E{(f),k) fc^sin0cos( 



+ 



A;' k'W^-k^ sin^ 



(93) 



which is the incomplete version of the expression in iDurandl fll953l ). The partial derivative with 
respect to the parameter m is: 



— 9„n(0,m2,A;) 
2m 



sm 



(1 — sin^ 0)2a/1 — k"^ 



(94) 



sm 



The integral (which can be obtained by an arduous calculation from the elementary decom- 
position of the denominator) is also given by Dieckmann (2011)1. We then have: 



m 



— k^ 



E{<P,k)-F{<P,k) + 



— k"^ 

— k'^ 



n(0,m^ A;) 



sin d) cos (h\/l — k"^ sin^ 



(m^ — A;^)(l — sin^ 



where m' = -\/l — is the parameter complementary of m. 



(95) 



D Partial derivatives with respect to the radius a 

The partial derivative of the first term in Eq. (|T5|) is 



da[5E{<p, k)] 



a + R 
6 



E(0, k) + 5dakxdkE{^,k) 



-(a + R+ ^(0, ^) - ^(^^ k). 



(96) 



where 5, and m are given by Eqs 
- 



and ([5j) respectively. For the second term, we have 



5 



-F{<f),k) 



^ ' ^ -F{^, k) + dak X 9fcF(0, k) 



{0" - R^){C + R^ - a") E{(t), k) 



1 



53 



2aC^ + (a + i?)' 



2a53 
2a 



k" 
Fi^,k) 



a' - R'){C + R' - a') F sin COS. 



2a53 



(97) 



sm 



The last term containing the elliptic integral of the third kind is not easy to determine, since 
both m and k depend on a. We have 



(9an(0, m?, k) = 5^11(0, m^, k) x Sq/c + 9^11(0, m^, /i:) x dam 



(98) 



^See [http : //pi . physik . uni-bonn . de/'dieckmain/l 
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and so 



a„n(0,m2,fc) 



E{(P,k) 



+ 



sin d) cos < 



+ 



m 



4 A;2 



— k'^ 



sin (/) cos — k"^ sin 



(m^ — /c^) (1 — sin^ 



— /c^ 
^(0,m^fc) 

dam 



dak (99) 



mm 



We also have 



dnk = - 



and 



where 



m^ — k'^ 

m 1 



m^ — A;2 



dam 



2am2 
52 1 i? + a 



_ {R - a)5^ 

" 2{R + a)am'^C ~ C^R-a 

e 



4 ; 2 2/1 /2\ T 2 2 

m — fc = m(l — m j — /e =m 



/2 



(52 



— m 



After some algebra, we find 

a + Rl 



da'n{(j),m ,k) 



a-RS^ 



3R-a + 



2Re 



n(0,m^ A;) 



a - R 1 r^/ , , X 

2 

m 

+ :: — sm COS ( 



a + Rf 

a + R R^ + 3a^ + C 
2{a-R)a {a - R^ + C 



(100) 

(101) 
(102) 

(103) 



Ei^,k) 



2aC 



k^ R^ + C-a'^ Jl - fc2 sin^ cj) a + R 
+ 



As we have 



da 



a-Rl 
a + R6 



k'^m^ ^1 _ A;2 sin^ (1 - m2 sin^ 0) a - i? 
2i?C^ 



1 
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3i?-a + 



;a + i?)2j ' 



(104) 



we deduce that the partial derivative of the third term in Eq. fllSp does not depend on 11. 
Actually, we get 



da 



a-RC 
a + R 5 



n(0,m^ A;) 



2a5 



2a(5 (a - i?)2 + (2 



(105) 



+ 



1 a — R sin cos 1 



+ 



2a5a + R ^1 - k^ sin^ m2/c 
m25 sin cos 1 — k'^ sin^ 



;,(i?2 + C2-a^ 



2a 



1 — m2 sin^ I 
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E Integration of the asymmetric term 

We first separate the term into two terms, namely 



a/1 — /c^ sin^ I 
f 1 — sin^ d 



We then get 
~2a 



■ sm cos ( 



yl — Psin^ 
'\ — rrfi sin^ d 



1 




k'^ sin^ 




1 




sin^ 




1 




sin^ 



1 - fc^ sin <; 
(1 — sin^ I 
1 

1 — sin^ d 



1 — sin^ I 



2i?5 sin d> cos < 



(106) 



fc^ rv? sin^ 
1 — sin^ i 

1 — ^ +^ 



(a + Vl-Psin^, 
1 



X 



— + 



1 — sin^ I 



+ 



2i? sin cos 
-4ai?sin2 

IRC' sin cos ( 



(107) 



[(a + i?)2 -4ai?sin2 0] v^^^ - 4ai? sin^ ' 

where 5^ — 4ai?sin^0 =PP'^. The integration of the first term (if we omit the quantity 
2R sin cos 0) with respect to a is: 



da 



la a/^^ — AaR sin^ 
Regarding the second term, we set 

a + R cos ! 



asinh 



a + R cos 20 



Ve + i^^si 



sm 



u 



a + Rcos 20 



As 



^/{a + R cos 20)2 + C + R'^ sin^ 



pp/ 



da 



1 - V(« + ^ cos 20)2 + C + R^ sin^ 20 PP 
the integral of this term now reads (if we omit the quantity 2('^R sin cos ( 



da 



a [(a + i?)2 -4ai?sin2 0]PP' 



du 



(2^2 + i?2 sin2 20 

1 



(108) 

(109) 
(110) 

fill) 



|Ci?sin20| 
1 



arctan 



c 



R sin 20 



a + R cos ! 
PP' 



-arctan 



C(a + i?cos 20) 
i?sin 20PP' 



Ci?sin20 

We can then write the second term in the right-hand-side of Eq. ffTSj) as a derivative, like 

a + R cos 



sin cos 0yT— Psir? 



2a 



'1 — m2 sin^ - 



dn < R sin 20 asinh 



+Carctan 



C(a + i?cos20) 
Rsin 20PP' 



112) 
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F Radial acceleration of the floating sector 

To determine the radial component of the acceleration g^°^^ = _ y^^fioat associated with Eq. 
we need the partial derivatives of k and m with respect to R. We have 



and 



ydRk 
k 



1 



2R5^ 
a — R 



113) 



(114) 



These relations can be deduced from Eqs.f ll6p and (ITTl) respectively by exchanging a and R. 
With the help of the Appendix [Dl we find 



k)] 



a + R 



6 



k) + 5dRk X dkE{^, k) 



;il5) 



a + R + 



2R 



and 



'R 



- R' 



F{<P,k) 



Ei<P,k) 



F{<P,k) 



e + a"- R' 
2R6 



Fi^,k), 



-2R6^ - {a^-R^){a + R) 
53 



- R^)iC + - R^) E{(f), k) 



a^-R' 



dRkxdkF{4>,k) (116) 



1 



2R6^ 
2RC + (a + Rf + 



k'' 

- R^){C + a' - R' 
2R 



Fi<P,k) 



- R^){C + - R^) A;^ sin cos. 



2R6^ 



fc'Vl - k'^ sin^ 



(117) 



where k' = \/l — k^, and 
'a-RC 



d 



'R 



a + R 6 



n(0, m, k) 



mr, a' + 3R' + e 



2/2(5' ' 2i?5 (a - /?)2 + (2 

1 a — R sin cos (/) fc^ 



:ii8) 



2R5a + R Jl-k^ sm^6 m^k'^ 



rri^S sin d) cos 1 — k^ sin 



2i? 1 — sin^ 

For the derivative of the two functions RH and (T, we respectively find 



dnRHi^) = H(4>) - 



sin 20 



C^cos20 - ai?sin^20 C^cos! 



+ R^ sin^ 



PP' 



and 



aRCT(0) = -c'sin; 



PP' 



X 



+ C^sin20 



i?2 sin^ 20PP'^ + (a + i? cos 20)2 

/ n A i?(i? + a cos 20) 
i? cos 20- (a + i? cos 20) ( 1 + ^ — ^7;;;2 ^ 

i? COS 20 



;ii9) 



(120) 



PP' 



iC + R^ sin^ 20)v/C^TR2 
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So the radial gradient of is 



float 



6 



R 



cos (t>\/l — k"^ sin^ < 



, , N I k \ „ , , , , m sin > 

z / 2 1 — m"^ sm i 



;i2i) 



i? sin 20 / cos 20 — ai? sin ^ 



cos 20 



(2 + i?2 sin^ 24 

-C^ sin 20 



PP' + 
pp/ 

i?2 sin^ 20PP'^ + C^{a + R cos 20)^ 



X 



i? cos 20 — (a + i? cos 20) , ^ , „ 
^ \ PP 



sin : 



/ + a cos 20) 

i? COS 20 



iC + R^ sin^ 2(f)) ^/WTC 



We see that this expression is fully compatible with the expression for daip'^^'^'^ given in the 
Appendix [B] for = f , since we get in this case 



float 



6_ 
R 



Eik) - ( 1 - y ) K(A:) 



(122) 



G Angular acceleration of the floating sector 

Regarding the angular acceleration dgip^"^^ of the floating sector, we have 



de = -dp, 



(123) 



meaning that we need in particular the partial derivatives of the elliptic integrals with respect 
to their amplitude /3. For any amplitude 0, we have 



a^F(0, k) 



1 



1 — k"^ sin^ I 



and 



9011(0, m, k) 



d^E{(j),k) = yl - Fsin^, 
1 



(124) 



(125) 



(126) 



[1 — m? sin^ 0) a/1 — fc^ sin^ 

None of the quantities k, 6, k and m depend on 0. The inverse hyperbolic term and the inverse 
trigonometric term in Eq.f lTS]) give respectively 

o I . . , a + i? cos 20 \ , . , a + -Rcos20 /^^^n 

Off) I sm20asmh — ^ _ ) = 2cos20asmh — ^ _ (127) 



y/C + R^ sin^ 



VC^ + sin^ 20 



o ■ 2o^^(C' + ^' + «^cos20) 
— 2 sm 20- 



PP'(C2 + R^ sin^ 20) 
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and 



Finally, we have 



C(a + i?cos2(/)) 
i?sin20PP' 



2CR PP' 



(a + - m2 sin^ 0)(C2 + gin2 

+ a cos 20 — (a + i? cos 20)/c 



2 sin^ cos^ I 



1 — k? sin^ I 



float 



PP' -R^ a-R 
—r- + 



2PP' a + /?2PP'(l -m2sin2 0) 

R{C + R^ + aR cos 20) R^/WTC 



(128) 



:i29) 



+ cotan 20 H{(f)) — sin^ 



PP'(C2 + i?2 sin^ 20) + R"^ sin^ : 

Rpp' r 

R + a cos 20 — (a + i? cos 20)/c 



iC + R^ sin2 20) t (a + ^)2(1 - sin^ ■ 



2 sm cos ( 
1 — A;2 sin^ , 



H Vertical acceleration of the floating sector 

We have dz = and so 

d^6E{(f), k) -- 
C(a2 - 



dc^^F{<p,k) 



6 



where k' = yl — W, and 
a-RC 



a + R 6 



n(0,m2,A;) 
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a-R C 
a + RS 



S ' 

' E{4>,k) 

k'^ k'^^l-k^ sin^ 



sin cos ( 



n(0,m2,A;) + — -2E(0,A;) 

A;^ sin cos 



(130) 
(131) 

(132) 



iv?k''^\/l — A;2 si 



sm 



The trigonometric term gives 



i?(PP'^ - C^)(a + i?cos20) 
sin^ 20PP'^ + C'^{a + R cos 20)^ 
C R{a + Rcos2(f)) 



9cCT(0) = T(0)-Csin20 



T(0) - sin; 



(133) 



PP' i?2 sin2 20 + (2 ' 



and the hyperbolic term gives 



d(;RH{(f)) 



R( sin : 



a + R cos 20 



R cos 20 



(2 + i?2 sin2 : 
Gathering terms, we finally find 



PP' + R^ sin^ 20 



float 



F{<P,k)-^U{<P,m',k) 
a + K 



+ T(0). 



(134) 



(135) 



This expression is also compatible with the expression for d(^ip'^^^'^ given in the Appendix [B] for 
0=1 since the term T(0) is — iCl^r, — |Clf or depending on the ratio R/a, according to 
Eq.dSH]). 
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I Potential of a rectangular sheet 



(Durand. 


1953; 


Chappell et al. 


2012) 



-iy'i 

+ (4 



y ) In 
y) In 
- x) In 
— x) In 





— T -\- \ / ( 






y) 


^ -j^ 


- X + a/ ix\ 


-X)2- 


Hy2 


-y? 


X2 


- X + a/ (4 


— x)2 - 


Hy'i 


-y? 




- X + a/ (x'l 


-X)2- 


Hy'i 


-yy 


y'2 


- 2/ + V (4 


— x)2 - 


Hy2 


-y? 


y'l 


- y + V (4 


— x)2 - 


Hy'i 


-y? 


y'2- 


- ?/ + v/(x'i - 


-X)2 4 


-{y'2- 


-y? 



y'i-y + V {x[ - x)2 + {y[ - y) 



(136) 



where r = xe^ + yey, and {x'i,y[), (xg > x[,y[), (x2,?/2 ^ y'l) and (x']^,?/2) denote the Cartesian 
coordinates of the four corners of the sheet. 
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